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Abstract 

Stolarsky [Proc. Amer. Math. Soc. 41 (1973), 575-582] showed a beautiful relation 
that balances the sums of distances of points on the unit sphere and their spherical cap L2- 
discrepancy to give the distance integral of the uniform measure on the sphere a potential- 
theoretical quantity (Bjorck [Ark. Mat. 3 (1956), 255-269]). Read differently it expresses the 
worst-case numerical integration error for functions from the unit ball in a certain Hilbert 
space setting in terms of the L2-discrepancy and vice versa (first author and Womersley 
[Preprint]). In this note we give a simple proof of the invariance principle using reproducing 
kernel Hilbert spaces. 



1 Introduction 

We consider the unit sphere 



4 + 



d+l 



embedded in the Euclidean space M'^^^, d>2. Let / : S'^ — )• C be a continuous function. Then 
we approximate the integral J^a fix) dad{x), where ad is the normalized Lebesgue surface area 
measure on S'^ (Jg^j d ct^ = 1), by an equal weight numerical integration rule 



1 



N-l 



fc=0 

where zq,...,Z]^^i E S'^ are the integration nodes on the sphere. In order to analyze the 
integration error committed by the approximation, we define a worst-case error by 



sup 

/eH,||/||<i 



f{x)dad{x) 



AT E f^^^) 



N 



k=0 



where Ti denotes a normed function space with norm || • ||. The rate of decay of the worst-case 
error depends on the function space and the integration nodes. For a fixed function space, the 
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worst-case error can serve as a quality criterion for different sets of integration nodes, meaning 
tliat the performance of a set of N quadrature points zq, . . . , ^at-i can be compared to another 
set of quadrature points by comparing the corresponding worst-case errors. Generahy, this 
only means that the integration error is smaller, but sometimes the worst-case error allows also 
a geometrical interpretation. 

Another quality criterion for points on the sphere exploits the potential energy, or more 
generally, the Riesz s-energy of configurations of points modeling unit charges which are thought 
to interact through a potential 1/|| • ||* (s / 0), where || • || denotes the Euclidean distance. (We 
refer the reader to the survey papers [10] and [13] and for universally optimal configurations to 
[_8j.) A particular instance is the (normalized) sum of distances (s = —1) 

^ X] II^A; - ^^ll> Zl, . . . ,ZAr e S'^. 

k,i=0 

It is well-known from potential theory (see Bjorck [5]) that this (discrete) sum of distances of 
optimal A^-point configurations approaches the associated (continuous) distance integral of the 
uniform measure on S'^ as — ?■ oo. In fact, any sequence of A^-point systems with this 
property turns out to be 'asymptotically uniformly distributed'; that is, the discrete probability 
measure obtained by placing equal charges at the points tends to the uniform measure in the 
weak-star sense. The difference 

N-l 



[ [ \\z - x\\dad{z)dad{x) - ^ y^\\zk 



k,e=o 

measuring the deviation between theoretical and empirical (— l)-energy quantifies the quality 
of points on the sphere (and, indirectly, their uniform distribution) using energy. It should be 
mentioned that the upper bound of correct order N^^^^^'^ for ([2]) (for optimal configurations) 
was obtained by Stolarksy [16] using his invariance principle and a result of Schmidt [H] on 
the discrepancy of spherical caps. The correct-order lower bound {N^^~^^'^) was established by 
Beck [3] using his Fourier transform technique. 

The spherical cap discrepancy measures the maximum deviation between theoretical and 
empirical distribution with respect to spherical caps as test sets. It can be used to compare 
point sets on the sphere with respect to their distribution properties. To introduce the concept 
of spherical cap discrepancy, we require some notation. A spherical cap centered at ic G S'^ with 
'height' t G [-1, 1] is the set 

C{x;t) = {z G S"' : {x,z) > t}. 
The family of all spherical caps is denoted by 

^ = {C{x;t) : a; G S^t G [-1,1]}- 
For a set J C R'^+i we define the indicator function 



1 if a; G J, 
otherwise. 



For a measurable set J C let 



cJd(J) = / lj{x) dad{x) 
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Using spherical caps we can define another quahty criterion for points on the sphere S*^ in terms 
of their distribution properties. One such criterion is the spherical cap L2-discrepancy, which is 
given by 

2 \ 1/2 



/ 1 N-1 

W-l 7§d 



dad{z)dt 



where P = {zq, . . . , zn~i}- 

We considered three, seemingly different, measures which can be applied to point sets on 
the sphere. It turns out that in some instances, the three measures are related to each other. 
Stolarksy's insight |16| was that the sum of distances of points on the sphere and the spherical 
cap discrepancy coincide. On the other hand, also the sum of distances of points on the sphere 
and the worst-case error coincide for a certain choice of function space, see [7] and also Sloan 
and Womersley [E] regarding a generalized discrepancy of Cui and Freeden [9] . 

In this paper we give a simple proof of these results based on reproducing kernel Hilbert 
spaces. We also provide some generalizations which follow from our approach. 



2 Reproducing kernel Hilbert space 

We define a reproducing kernel Hilbert space using the general approach of [12^ Ch. 9.6]. 
For x,y eE,'^ we define the function K^^ : S'' x S"^ ^> M by 



Kc^{x,y)= I / lc(z;t){x)lciz;t)iy) dcTdiz) dt. 
-1 J§<^ 



(3) 



Since lc{z;t)(x) = lc(x;t){z), we also have 

Kv{x,y)= [ [ lc(a;-t){z)lc{y;t){z) dad{z) dt. 
J -I Js<i 

The function is obviously symmetric, i.e. we have K<^(x,y) = K^^(y,x). Further, let 



, aN-i G C and xq, . . . , a;jv_i G S"'. Then we have 



N-l „i „ N~l 

^ akaiK^^{xk,xi) = / ^ X] «fcO?lc(a:;i)(aJfc)lc(a;;t)(a^^) do-d(a^) dt 



k,e=o 



k/=0 
N-l 



k=0 



dad{x) dt 



> 0. 



Thus, the function K^^ is symmetric and positive definite. By this implies that K^^ is a 
reproducing kernel. It is also shown in [2j that a reproducing kernel uniquely defines a Hilbert 
space of functions with a certain inner product. Let H'^ = 'H{K<^,E'^) denote the corresponding 
reproducing kernel Hilbert space of functions / : S"^ — )• M with reproducing kernel K<^. 

We consider now functions /i , /2 : S'^ — )• C which permit a certain integral representation. 
Let c/1,52 : S'' x [-1,1] C with 51,52 G IL2(S'^ x [-1,1]) and 



[ [ 9i{z;t)lc{z;t)ix) daaiz) dt, 



i = 1,2. 



(4) 
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Notice that for any fixed y G S'^ the function K^g{-, y) also is of this form, where the function g 
is given by lc(z;t)(y) (considered as a function of z and t and where y is fixed). For functions 
of this form we can define an inner product by 

{fi,f2)K^=[ [ gi{z;t)g2{z;t)dad{z)dt. (5) 

Let y €z S'^ he fixed. With this definition we obtain 



{fi,K^^{-,y))K^ 



I I 9i{z;t)lc(^.t){y) dad{z) dt = fi{y). 



By [2], the inner product in is unique. Therefore, functions /j, which are given by @ and 
for which (/j, fi)K^ < cO; are in and ([5]) is an inner product for those functions in H^g. 
Consider now the reproducing kernel Keg. We have 

/I f-vam{{x,z),(y,z)} 
^'^C(z;t){x)lc{z;t){y) = J ^ dt = l + min{(a;,2;),(y,2;)}. 



Thus 



We have 



and 



Kv{x,y)= [ [ lc(z;t){x)lciz;t){y) dad{z) dt 
= 1+1 min{{x,z),{y,z)] dcFd{z). 



mm{{x, z),{y,z)} = ^ [{x,z) + {y,z) -\{x-y,z) 



/ {x,z) dadiz) = 0. 
X = y, then we therefore obtain 

/ mm{{x, z),{y,z)} da d{z) =0. 

Let now x ^ y. Then 

/ mm{{x,z),{y,z)} dad{z) = -I- I \{x - y, z)\ dad{ 



„ 1 

\x — y|| - 

2 J§d 



x-y 

\x - y|l 



dad{z). 



The last integral does not depend on the unit vector {x — y)/||a; — y|| by rotational symmetry. 
Thus, we have (cf. Appendix |A|) 

2 7§d a ijJd fid\p'^) a V^r(d/2) ^/2^^d 

(Here, ujd is the surface area of S'^, B'^ is the unit ball in M°', T-Ld is the d-dimensional HausdorfF 
measure normalized such that the d-dimensional unit cube [0, 1)'^ has measure one and T{z) is 
the Gamma function.) 

Therefore we obtain the following closed form representation: 

Keg{x,y) = l-Cd\\x-y\\. (7) 

The reproducing kernel has the following properties: for y, z G S*^ we have 
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. l-2Cd<K{x,y)<l- 

• K{x, y) = I — 2Cd if and only if a; = —y (clearly, < 2Cd < 1); 

• K{x, y) = 1 if and only if a; = y; 

Note that the Karhunen-Loevy expansion of the function — y|| is based on ultraspherical 
harmonics. Hence the eigenfunctions of are the ultraspherical harmonics. The corresponding 
eigenvalues are also known. Therefore, the functions in Tii^ can be expanded using ultraspherical 
harmonics and the inner product can be written using the coefficients of such an expansion. See 
[7] for these results. 



3 Worst-case error 



Let ll/llii'^ = V (/, f)K'^ denote the norm in W^^. Then we define the worst-case error for a 
quadrature rule Qn given in ([T]) by 



eiW^, Qn) = sup 



/ f{x)dx-QN{f) ■.f€n^^,\\f\\K^<l]. 



Let / G Ti'^. Then, by the reproducing kernel property f{y) = {f,K<^{-,y))K^ for y G S"^, 
and, since the integration functional / i— J^^ f dad is bounded on "H"^ and has the 'representer' 
Jgd Kc^(-, z) dad{z), one can write 

N-l 

/ fix) dadix) - ^ E •^(^'^) = {f,nn^^,QN; ■))k^ , 

where the 'representer' of the error of numerical integration for the rule Qn for functions in Ti^^ 
is given by 



1 

'JI{W^,Qn;x) = / K^^{x,y) dadiy) - — K<^ix,x 
•^^ fc=0 



X £ 8". 



The Cauchy-Schwarz inequality yields 

1 

/ /(a;) dCTd(a;) - — V /(a;; 



k=0 



{f,n{nv,QN;-))K^ < \\f\\K^mm'^,QN;-)\\K^. 



In particular, equality is assumed in the last relation when taking / to be the 'representer' 
TZ{H^, Qn] •) itself. It follows that 



e{'H^^^,QN) = m'H^^,QN;-)\\K^ 



1 

/ K.4.,x) dadix) --Y,K^d- 



(8) 



Expanding the square of the worst-case error and substituting the closed form of the reproducing 
kernel we arrive at the well-known representation 

[e{n<^,QN)f = {n{W^,QN;-),n{W^,QN;-)) 

N—1 /V— 1 

= / / K<g{x,y) dad{x) dadiy) - tt / K^^{x,Xk) dad{x) + ^ V K^^{xk,xe) 
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\x - y\\ dad{x) dad{y) 



1 



N-1 



k,l=Q 



xe\ 



(9) 



This shows how the square worst-case error in our reproducing kernel Hilbert space is related 
to ([2]). In the next section we show how the worst-case error e(7i^^, Qn) in our reproducing kernel 
Hilbert space relates to the L2 spherical cap discrepancy. 



4 Spherical cap discrepancy and Stolarsky's invariance principle 

Using the integral representation of the reproducing kernel ([3]) we have 



Kc^{x,y) dad{y) 



'^C{z-t){x)ad{C{z]t)) dad{z) dt 



and 



^ N-l r-l r. ^ N-l 

^ Kcg{x,Xk) = I,nY1 ^C{a^;t)iz)lc{a^,;t){z) dcTd(z) dt. 

Thus, the 'representer' of the error of numerical integration is of the form (jl]); that is 







n{'H'ff,QN;x) = [ 1 


lc(z;t)(a3) 


J -I J 



N-l 



k=0 



dadiz) dt. 



Therefore, using the inner product representation (l5|) in ([8]), we obtain 

2 \ 1/2 



/ 1 N-l 

e{W,,QN)= / / ad{C{z;t))--Y^ 



C{z;t)\Xkj 



dad{z) dt 



(10) 



stating that the worst-case error of the numerical integration formula Qm in dH) in the considered 
Sobolev space setting equals the so-called spherical cap L2-discrepancy of the integration nodes. 

Combining ([9]) and (jlOp . we arrive at Stolarsky's invariance principle for the Euclidean 
distance on spheres. 

Proposition 1 (Stolarsky [16j ) Let xq, . . . ,xj\[^i G S'^ be an arbitrary N point configuration 
on the sphere S'^. Then we have 

2 

dcJd(z) dt 

k=0 



N-l 



N-l 
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1 \ f'^ r 1 

= / / 11^ - y\\ do-d(a3) d(Td{y). 

The L2-discrepancy of an A^-point configuration on S'^ decreases as its sum of distances 
increases and vice versa. The right-hand side is the distance integral of the uniform measure ad 
on the sphere S*^ which is the unique extremal measure (also known as the equilibrium measure) 
maximizing the distance integral 



:= / / \\x - y\\ d ii{x) d ii{y) 



over the family of (Borel) probability measures fj, supported on S'^. For the potential theory of 
the generalized distance integral we refer to Bjorck 
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5 A weighted reproducing kernel 

The above results can be generalized by introducing a weight function. Let v : [—1,1] — ?• M 
satisfy v{t) > for all t and which has an antiderivative, which we denote by V. Then we define 
the reproducing kernel with weight function v as follows 



Kgf^y{x,y)=[ v{t) [ lc{z;t){x)lc{z;t){y) dad{z) dt, x,yeS'^ 

J-l J§d 

For functions represented by integrals 

fi{x)=[ [ gi{z;t)lc{z;t){x) <i(^d{z) dt, i = 1,2, 
the corresponding inner product is now given by 



(12) 



1 



v{t) 



gi{z;t)g2{z]t) dad{z) dt. 



The reproducing kernel can be written as 

K.g^^{x,y)= f V{in\n{{x,z),{y,z)]) dad{z)-V{-l), x,y£§'' 



(13) 



For certain weight functions v, this expression may have a concise form. This reproducing kernel 
defines a reproducing kernel Hilbert space 

The 'representer' of the error of numerical integration for the rule Qn for functions in H^g^^ 
takes on the form 



1 Jn- 



lc{z;t){x)v{t) 



ad{C{z;t)) -j^Yj ^C{x,;t){^ 



k=0 



dad{z) dt. 



We claim that Kcgy[x,y) is a function of the inner product {x,y), cf. Appendix [Bj Using the 
same approach as before we obtain 



N-l 



[e{'H'g^y,QN)]^ = ^ K'^,v{^k,xe) - [ / K<g^^{x,y) dad{x) dadiy). (14) 
This worst case error can also be expressed in terms of a weighted discrepancy measure: 

N-l 



ei'H^g,v,QN) 



■^^ fc=0 

1 ^ N-l 

J /^^^ Jsd ^'^^'-^^^'^^^ ~ iV ^ lc(a:;t)(aJfc) 



= \\'^{T~(-'^,v,Qn; ■)] 

2 \ 1/2 

dad{x) dt 



Using p4p and (|15p we obtain the weighted version of Stolarsky invariance principle. 



(15) 



Theorem 1 Let Xq, . . . ,X]^-i £ S'^ be an arbitrary N point configuration on the sphere S'*. Let 
K^g^v be the weighted reproducing kernel given by . Then we have 



N-l 



N-l 



^ K^g,v{xk,xe) + / v{t) I ad{C{x;t)) - '^C{x;t){xk) 
k,e=o '— " 

K^g,v{^,y)dad{x)dad{y). 



k=0 



dad{x) dt 
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The double integral above can be expressed in terms of the weight function, see (|20p . 

In [5], Stolarsky's (general) invariance principle is extended and used to get bounds for the 
spherical cap discrepancy, see also [6] . Stolarksy [T7] also extended his principle to certain metric 
spaces arising from measures. 

Stolarsky ^16j introduced the function 



p{x,y) 



max{(a;,y),(y,t)} 



§<* Jvam{{x,y),{y,t)} 



g{u)dudad{z), x,ye8 



which becomes a metric if the kernel g (integrable on [0, 1]) is positive but the proof of the 
corresponding invariance principle 



N-l 



p{x,y) dad{x) dad{y) 



k,e=o 



N-l 



k=0 



dad{x) dt 



(16) 



making use of Haar integrals over the special orthogonal group SO{d + 1) does not require it. 
Note that for 5=1 the function p{x,y) is a constant multiple of the Euclidean distance. With 
some care one may even consider g{x) = 1/(1 — x^). 

It is well-known that a reproducing kernel K{x, y) induces a distance (metric) by means of 



d{x, y) = dxix, y) = \J K{x, x) - 2K{x, y) + K{y, y). 
For example, the reproducing kernel d?]) yields 

dK^^.y) = \/c'd\/\\x - y\\. 



In general, for the symmetric weighted kernel K^^^v{x,y), which does only depend on the inner 
product {x,y), it follows that (a G S'^ fixed) 

If 2 1 1 2 

By Theorem [1] on arrives at 



N-l 1 N-l 

[d^^,.(^, y)] ' + 2 / v{t) / ^ ad{C{x; t)) - _ ^ 1 
kl=0 •'"^ •'^'^ k=0 



^C(x;t)[Xkj 



dad{x) dt 



which should be compared with (I16p . 



[dr^A^^y)]'^ dcrd{x)dad{y), 



Acknowledgement: The first author is grateful to the School of Mathematics and Statistics 
at UNSW for their support. 

A Auxiliary results 

The normalized surface area measure ad on S'^ admits the following decomposition 

,2\d/2-l 



dadiy) = ^{l-t'y 

^d 



dtdad-i{y*), y = {Vl-t^y*,t)eS'' 



(17) 
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where t E [— 1, 1], y* G S'^ ^ and ujd denotes the surface area of S'' (cf Miiller [II]). (By definition 
iUiP) = ^) where p is the North Pole of S°'.) Thus, by rotational symmetry, the integral of a 
zonal function f{{z, •)), z G S'^ fixed, with respect to ad reduces to 

/ fi{z,y))dad{y)= [ f{{p,y))dad{y) = '^ f f{t){l-fY'^-'dt. 

J^d J§d J-1 

Proof. [Proof of relations ^] One gets 

Cd = U \ip,y)\aa4y) = l^ f \t\ {l - t^f^-Ut = [\l - t^^'-' 2Mt 

2 Jgd 2 J -I 2 i^d Jo 

_ lood.i _ i r(((i+l)/2) _ HdiM'^) 1 
The second equality follows from 

l = ad{§^) = ^ f\l_t'f'-Ut = ^ [\'/'-'{l-vf'-'dv = ^B{l/2,d/2), 

^d J-l ^d Jo '^d 

where B(a,6) = r(a)r(6)/r(a + 6) is the beta function. The third equality follows from the 
well-known formulas for the volume of the unit ball in M*^ and the surface area of S*^. The 
asymptotics follows from the asymptotic expansion of a ratio of Gamma functions (cf. [T]). □ 



B The weighted reproducing kernel 

Next, we investigate the weighted reproducing kernel ()13p in more detail. In particular, it will 
be shown that the kernel Kc^^^{x,y) is a function of the inner product {x,y). 

On observing that {x, z) < {y, z) if and only if {y — x,z) > we may write for x ^ y 

K^^,v{x,y) = A^^y{x,y) +A^^y{y,x) - V{-1), 

which immediately shows symmetry of the reproducing kernel, where 

A^,,{x,y)= [ V{{x,z))l[o^^]{( y ~ z))dad{z). (18) 

Js'' \ \\y - x\\ I 

By abuse of notation we set (note that (a;, y) = u) 

z = tx + Vl-t^ z*, -1 < t < l,z* e S'^-\ 
y = ux + y/l-u'^y*, -1 < t < l,y* £ S'^'^ . 

In this way x will be the 'North Pole' in the decomposition (117p and we obtain 

A,A^,y) = ^ f v{t)\ [ l[o,i](//^,z\w_i(z*)| {l-t^'-'dt. 

^d J -I Us<i-i \\\y / J ' 

The indicator function in the inner integral is a zonal function depending on w = {y*, z*) only. 
Thus, we apply again (I17p with y* as 'North Pole'. That is 



9 



where the inner product evaluates as 

(^..) = v/T3F/I±E«I..-yiIE5[i. s^.^^. (19) 

Proceeding similarly for A^^^v{y,x), one sees that, indeed, A'^^v{y,x) = A^^v{x,y). Further- 
more, A'^^y{x,y) depends only on the inner product {x,y) which in turn implies that the 
reproducing kernel K<^^y{x,y) is a function of the inner product {x,y). 

The right-hand side in (llOp describes a line which stays strictly between the levels —1 and 1 
for V in [—1,1] by the left-hand side in (jl9p . Further analysis gives that the indicator functions 
is o ne (i) i f t < -y/{l + u)/2 and -1 < < 1, or, (h) if -y^(l + n)/2 < t < y'{l + u)/2 and 
{t/Vl -t'^)^/{l - u)/{l + u) <v<l, and zero otherwise. This leads to 



J -I 



+ ^ f^"^^ I(i-.w)/2((d - l)/2, {d - l)/2) (1 - e f^'Ut, 
where, when using u = {x, y) = cos (0 < (/> < vr) and t = cos V', one has 



-cos{4>/2), J,^=tan((^/2), - ' 



2 S l + u " Vl + nVr~t2 cot((?i/2)' 

The change of variable ^ = x{t) yields 

V(l+«)/2 



^u^/i-ny/^ r^^^^^ de 

Wrf VI + J-1 Jl^^f2 



S + e^ (Its 

where we make use of the regularized incomplete beta function 

h{a,b) =B^{a,b)/B{a,b), B^{a,b) = [ t^"^{l - tf-^ dt, a,b>0. 

Jo 

We compute the following integral (using (fT8|) ): 

/ Ar^^^{x,y) dadiy) = V{{x,z)) l[o,i]((||^ — ^,z)) da^iy) da^iz). 
Jsd Jsd Jsd \ \\y - x\\ I 

The inner integral is one if z is in the half-sphere centered at — and zero otherwise. Hence, 

/ A^,,(^,^)da,(^) = ^ l\^-t){\-e)'"^-^dt. 

J§d ^d Jo 

Since g'{t) = {u;d-i/oJd){l - t^fl^-^ for g{t) = (1/2)1^2(1/2,^/2), < t < 1, integration by 
parts gives 

/ A^^,{x,y)dad{y) = \v{-l) + \ f v{-t)lt2{l/2,d/2)dt. 
Js'' ^ ^ Jo 



It follows that 




K.g,,{x,y)dad{x)dad{y) = 2'^ [\{-t) {l - t^f^ 'dt-V{-l) (20a) 

= / v(-t)It2(l/2,(i/2)dt. (20b) 
Jo 
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